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Low-energy kink in the nodal dispersion of copper-oxide superconductors: Insights 

from Dynamical Mean Field Theory 
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Motivated by the observation in copper-oxide high-temperature superconductors, we investigate 
the appearance of kinks in the electronic dispersion due to coupling to phonons for a system with 
strong electronic repulsion. We study a Hubbard model supplemented by an electron-phonon cou- 
pling of Holstein type within Dynamical Mean Field Theory (DMFT) utilizing Numerical Renormal- 
ization Group as impurity solver. Paramagnetic DMFT solutions in the presence of large repulsion 
show a kink only for large values of the electron-phonon coupling A or large doping and, contrary 
to the conventional electron-phonon theory, the position of such a kink can be shifted to energies 
larger than the renormalized phonon frequency uiq. When including antiferromagnetic correlations 
we find a stronger effect of the electron-phonon interaction on the electronic dispersion due to a 
cooperative effect and a visible kink at wj, even for smaller A. Our results provide a scenario of a 
kink position increasing with doping, which can be related to recent photoemission experiments on 
Bi-based cuprates. 

PACS numbers: 74.72.-h,71.38.Mx,71.10.Fd,71.20.Tx,71.10.-w,63.20.Kr 



I. INTRODUCTION 

One of the manifestations of a strong coupling between 
condensed matter electrons and lattice vibrations is the 
occurrence of kinks (abrupt changes of the slope) in the 
electronic dispersion relation. Such kinks occur at a cer- 
tain energy associated with a phonon mode and their 
strength is directly related to the electron-phonon cou- 
pling constant A. Therefore the observation of a kink, for 
instance, in angle resolved photoemission (ARPES) data, 
is often taken as an indication that a system possesses a 
strong electron-phonon coupling and its properties are 
likely to be strongly influenced by this. However, also 
coupling to other bosonic modes possibly of purely elec- 
tronic origin, can have the same effect on the dispersion 
relation. Moreover, in addition to an electron-phonon 
coupling different interactions important at low energy 
can compete or cooperate with it, which may modify the 
picture. Therefore, the experimental observation must 
be analyzed carefully before conclusions about the origin 
of a kink can be made. 

Prominent examples for the observation of kinks come 
from the field of cuprate superconductors. Studies by 
Lanzara et alX collected evidence for this in different 
compounds. One interpretation of these experiments 
is that the dispersion kink along the nodal direction 
of the Brillouin zone is associated to a longitudinal in- 
plane bond-stretching mode in which two of the four oxy- 
gen atoms surrounding the copper move inwards ("half- 
breathing" mode)^.. The same phonon mode indeed dis- 
plays a sizable renormalization in inelastic neutron scat- 
tering measurements 3 , pointing towards a strong cou- 
pling to carriers. This conclusion is further supported 
by Iwasawa et al£, who carefully studied the effect 
of an 16 0- 18 isotope substitution and concluded that 
their observation is compatible with the "half-breathing" 



mode. On the other hand, recent work by Dahm et 
al£- emphasized the connection between the nodal kink 
and spin-fluctuations, an alternative scenario £~— Hence, 
presently no general consensus exists on this important 
issue of high-temperature superconductors. 10 A review 
article covering numerous aspects of this topic has been 
recently published by A. S. Mishchenkc 11 . 

One of the aspects of the nodal kink which is particu- 
larly puzzling is the doping dependence of its position. 
Early studies indicated that the nodal kink occurs at 
a more or less constant energy position 12 . A more re- 
cent study by Kordyuk et al. 13 showed however that in 
LSCO the kink position u>k first increases with doping, 
exhibits a maximum around optimal doping, and eventu- 
ally decreases again, in a dome- like shape. In Bi-2212 the 
situation is even more intricate as the non-monotonic de- 
pendence is replaced, at low temperatures, by a peculiar 
increase upon increasing doping. Such a pronounced dop- 
ing and temperature dependence is a very strong indica- 
tion that the nodal kink in cuprates cannot be described 
within the textbook picture of conventional metals with 
electron-phonon interaction. 

The standard description of effects of electron-phonon 
coupling is based on perturbation theory in the electron- 
phonon coupling constant assuming that electrons can be 
described as a weakly interacting Fermi liquid 14 . Then, 
to lowest order in the coupling the magnitude of a kink, 
TOfc, in the electronic dispersion relates to the dimen- 
sionless electron-phonon coupling A as m^ = 1 + A. A 
clear kink, however, points to a relatively strong electron- 
phonon coupling, for which one expects effects beyond 
the lowest order diagrams and possibly vertex corrections 
to become very important. Vertex corrections are beyond 
the Migdal-Eliashberg formalism, i.e. beyond the most 
common diagrammatic theory for the electron-phonon 
interaction* 1 ^—. If on top of that one takes also the 



electron-electron interaction and Mott physics into ac- 
count, as in any minimal model for high-temperature 
cuprate superconductors, a quantitative description of 
the nodal kink within conventional theories seems incon- 
ceivable. 

Substantial progress in understanding strongly cou- 
pled electron-phonon systems has been developed in 
recent years based on dynamical mean field (DMFT) 
methods 1 ^ and cluster extensions! 9 . The advantage of 
these approaches is that electron-electron and electron- 
phonon interactions are treated non-perturbatively and 
on equal footing 2 ^— . One major insight of these studies 
is that the electron-phonon interaction in strongly cor- 
related metals cannot be understood with conventional 
tools. In particular, the polaronic signatures in pho- 
toemission spectra are more complicated than what sug- 
gested by perturbation theory on top of a Fermi liquid 
solution 14 . This is also in line with the conclusions of 
other studies which do not rely on the dynamical mean 
field approximation 2 ^—. 

Even though DMFT misses a crucial part of the 
non-local physics of the two-dimensional copper-oxygen 
planes, it is interesting to see whether some aspects of 
the nodal kink puzzle posed by recent experiments can 
be understood already within the purely local picture 
of the interplay between the electron-electron repulsion 
and the retarded phonon-mediated attraction given by 
DMFT. This is the motivation for the present study, 
which focuses on a description of the kinks in the elec- 
tronic dispersion due to phonons in the presence of strong 
electronic correlations. We employ the numerical renor- 
malization group (NRG) 31 ' 32 as DMFT impurity solver. 
NRG is known to have a good low energy resolution at 
low temperatures 3 ^—, and one can calculate spectral 
functions directly on the real axis, which makes it the 
method of choice for the analysis of low-energy kinks. 

We first analyze paramagnetic (PM) solutions in the 
vicinity of half-filling and zero temperature. In order 
to see the effect of antiferromagnetic (AF) correlations, 
which are expected to be important in the cuprates, we 
also perform calculations in the commensurate AF state. 
The latter have been shown to successfully mimic the ef- 
fect of antiferromagnetic correlations and exchange cou- 
pling Jj£r21 Away from half filling the Hubbard model 
can show phase separation or incommensurate order, 
hence, the analysis of the results must be done very care- 
fully. 

Our study is based on a model Hamiltonian, Hubbard 
Holstein model, of the form 
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c\ a creates an electron at lattice site i with spin a, and b\ 

a phonon with oscillator frequency wo , ftj,cr = Ci a Ci^. We 
include nearest neighbor hopping tij with magnitude t. 




G 



FIG. 1: Lowest order diagram for the electronic self-energy 
due to electron phonon coupling. 



The electrons interact locally with strength U, and their 
density is coupled to a local phonon mode with coupling 
constant g. We have set the ionic mass to M = 1 in Eq. 
([1]). We define A = po2g 2 /u!Q in terms of bare parameters 
where po is the electronic density of states at the Fermi 
energy. 

For the no n- interacting density of states po(e) we take 
the semi-elliptic form po(e) = 2\/D 2 — e 2 /(ttD) 2 . In all 
the results we present here we take the value W = 2D = 
At = 4 for the bandwidth. The DMFT and NRG cal- 
culations are carried out as described in detail in Refs. 

The Paper is organized as follows: In section [H] we 
discuss the occurrence of kinks in the electronic disper- 
sion in the normal, paramagnetic state. We first look 
at results for U = and then analyze the situation for 
large U. In section Hill we show results when the effect of 
antiferromagnetic correlations are included in the calcu- 
lations. In Section HVl we conclude and relate our results 
to the experimental observations. 



II. KINKS IN THE NORMAL STATE 

Before analyzing the strong coupling situation we 
briefly recall the description of electron-phonon kinks in 
the weak coupling theory (for a review see, e.g., Refs. 
13811391 ). In the paramagnetic state the Green's function 
reads, 



G fc (w) 



1 



£k + M _ s ( w ) ' 



(2) 



where we assume a fe-independent self-energy. 4 - Gener- 
ally, the self-energy has contributions due to electron- 
phonon coupling g and electron-electron repulsion. 
The electron spectral function is given by pk{u) — 
— ImGk(uj)/ir and the full (interacting) dispersion rela- 
tion E k = E(ek) by the maxima of pk{oj). This corre- 
sponds usually to the poles of Eq. ([5]), i.e. the solution 
of 



E k -£ k + n-ReE(E k ) = 0. 



(3) 



We define a kink as (more or less) abrupt change in the 

slope of E{ek), which consequently must be found in 

ReS(w). 

In a calculation to order g 2 (depicted in Fig. QJ one finds 

that the Re£(u;) is roughly linear for small |w| <C Wo> 



ReS 9 (w) = —Aw, strongly peaked at lj ~ ujq (logarith- 
mic divergence), and small for |w| 3> u>q. This yields for 
the low energy dispersion E l k = £fe/(l + A) for \Ek\ <C w 
and E^ — e k for |i?fe| ^$> ujq. The magnitude of the kink is 
characterized by the ratio of the high energy over low en- 
ergy slope, which gives rrik = 1 + A as mentioned earlier. 
This quantity is equal to the inverse of the quasiparticle 
renormalization factor Z and for a fc-independent self- 
energy to the effective mass m* /mo = Z~ x . The free 
phonon propagator is sharply peaked at \uj\ = lo . Due 
to the coupling to the electrons the peak can be shifted to 
a renormalized value Wq and broadened. For the phonon 
properties, we consider the function B(u) = ((b;w)) u , 
which can be calculated in the NRG from the matrix ele- 
ments and excitations. The corresponding spectral func- 
tion is pb{u) — — ImB(uj)/ir. 



A. Results at U = 

We first want to compare our DMFT-NRG results with 
the simple minded theory and consider the case U = 0. 
We use luq — 0.2t, and a filling factor of n — 0.9. In Fig. [5] 
(a) an example of the Sk resolved spectral function /3fe(w) 
is plotted for A = 0.29. 




kinks in the dispersion occurs (see Fig. H](b)). At this en- 
ergy the imaginary part of the electronic self-energy sets 
in, and consequently the electronic spectrum becomes 
much broader. 

The dispersion relation E^ is obtained numerically 
from the peak of position of pk{u) as a function of e k . It 
is shown for two values of A in Fig. [3] (b). We find the 
picture of a kink as inferred from the lowest order g 2 dia- 
gram, i.e., the low energy dispersion is renormalized and 
at the (renormalized) phonon scale wj it reverts back to 
the non-interaction slope. The magnitude of the kink mu 
can be inferred from the low-energy slope. For A = 0.29, 
we find the value m^ = 1.75, whilst the corresponding 
value from the second order theory is with ink = 1-29 
substantially smaller. Comparing with part (a), we can 
also see how the discontinuity in the dispersion relation 
comes about through the change in maximum position. 

The corresponding behavior of the real part of the self- 
energy, S R (w), is shown in Fig. [2] (c). We can see how 
the low energy slope increases with A together with a 
peak at the renormalized phonon frequency u!q, smaller 
than the bare value luq = 0.2. The renormalized phonon 
frequency ujq as a function of A is shown explicitly in Fig. 
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FIG. 3: (Color online) PM calculation: The renormalized 
phonon frequency ujq as a function of A for U — and U = 6. 
The inset shows the phonon spectral function pbiy) for [7 = 6 
and A = and A = 1.34. 



FIG. 2: (Color online) PM calculation: (a) The Eh resolved 
spectral function for U — and A = 0.29. (b) The dispersion 
relation E^ as a function of £*. for U = and different values 
of electron-phonon coupling. The dashed line is for free elec- 
trons, Eh — £fc. (c) The real part of the self-energy S r (oj) 
(with E R (0) subtracted) as a function of to for U — and 
various values of A. (d) The DMFT result for the effective 
mass (circles) compared with the second order perturbation 
theory result 1+A (straight line). 



We can see sharp quasiparticle peaks at low energy and 



broad peaks from certain energy on, uj k 



where the 



A comparison for 1 + A with the DMFT results for a 
series of values can be found in Fig. [5] (d). There are 
quantitative deviation already at relatively small values 
of A and very notable ones for larger A J^ Hence, for neg- 
ligible U a very large kink is obtained already for values 
of the order A ~ 0.4. For larger values of A a bipola- 
ronic instability occurs where local pairs are formed. At 
half filling above a critical coupling stable solutions of a 
bipolaronic insulator can be found as in previous DMFT 
calculationSi 16 i 17 i 20 i 41 However, in the doped system the 
calculations become unstable, and results oscillate be- 
tween solutions with small and large double occupancy. 



Paramagnetic DMFT calculations can not provide a suit- 
able ground state solution in these situations. 



B. Results for finite U 

We now turn to the situation with finite Coulomb re- 
pulsion, where we have contributions to S(w) from both 
interaction terms. Therefore, already for A = electrons 
are renormalized due to scattering processes, which leads 
to an effective mass enhancement and finite lifetime. For 
small values of the Coulomb interaction, U < D, qualita- 
tively similar results to the ones discussed for U — can 
be found for the dispersion relation and the occurrence 
of kinks. 

For larger values of U the situation changes. One 
important thing to consider is that we have Holstein 
phonons, i.e. lattice vibrations which couple to local fluc- 
tuations of the electronic density. Such phonon modes 
are in direct competition with the Hubbard repulsion, 
which does exactly the opposite: it tries to freeze density 
fluctuations out. As a consequence, larger values of A 
(compared to U — 0) are necessary for observing effects 
on the dispersion. Furthermore, for values of U rele- 
vant for strongly correlated materials, U > W, the high 
energy slope does not revert to the non-interacting one 
for £fe € (—D,D). Thus, when characterizing the kink it 
does not make sense to compare with the non-interacting 
slope. Instead one can compare the high and low energy 
slope within one calculation. 

In order to be in a parameter regime roughly suitable 
for cuprate physics, we choose n=0.9 (10% hole doping) 
and U=Q > U C 2, U c i — 5.88 being the critical value for 
the metal-insulator transition at half-filling (n=l). 42 We 
first report the local electronic spectral function for A = 
and A = 0.86 in Fig. 5] The well-known three-peak 
structure can clearly be distinguished, with lower and 
upper Hubbard bands and a quasiparticle band at the 
Fermi level cu—0. The spectrum is not symmetric because 
of the finite hole doping. 

One can observe in Fig. U that the central quasiparti- 
cle peak is relatively little affected by the presence of 
the electron-phonon coupling. In fact, the effect of the 
electron-phonon interaction on the low-energy quasipar- 
ticle appear just as a small "de-renormalization"; in other 
words a decrease in the effective mass m* (see Fig. [S}i for 
small values of A) . This result has been already observed 
at half-filling within DMFT using Lanczos as impurity 
solver 21 . It can be understood from the physics of the 
effective impurity model and the Kondo effect connected 
with the quasiparticle peak. It comes from the fact that 
the electron-phonon interaction effectively reduces the 
bare repulsion of the system, making the Kondo cou- 
pling somewhat larger. Such a physical mechanism can 
also be described within perturbation theory, assuming 
that the scale U associated to the separation between the 
Kondo peak and the Hubbard bands is the largest, com- 
pared to t and o;q 43 . It is remarkable that we observe this 
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FIG. 4: (Color online) PM calculation: The local electronic 
spectral function for U — 6 and n = 0.9. 



de-renormalization effect so clearly even away from half- 
filling, i.e. where the abovementioned separation of en- 
ergy scales is no longer clearcut (at least for the removal 
part of the spectrum). Moreover this result convincingly 
shows how this kind of "protection" of the low-frequency 
part of the spectrum from phonon effects characteristic of 
the paramagnetic DMFT solution is present also within 
the NRG solution. At higher energies instead, one can 
see differences between the solution with and without A. 
Here, contrary to the low-frequency regime, the resolu- 
tion of our NRG calculation is not accurate enough to 
resolve the details of the spectra, because of the conven- 
tional high-energy broadening used. 

Let us now turn to the discussion of the dispersion 
kink. In Fig. [5] we show a collection of results for the 
situation with the same U and n and different values of 
A. The low energy part of /9fe(w) for A = 0.86 is shown in 
Fig. [5] (a) and the corresponding dispersion relation in 
panel (b). 

Comparing the result for A = and A = 0.86, we can see 
how the whole behavior is modified due to the electron- 
phonon coupling even though no clear kink is visible. 
We just get a continuous change of the slope for rela- 
tively large values of the coupling. We can understand 
this in light of the arguments given above: the Holstein 
electron-phonon interaction can not easily induce charge 
fluctuations due to the direct competition with the Hub- 
bard U . Even though the finite doping reduces the in- 
fluence of the Mott physics compared to half- fillin g 20 ' 44 , 
the whole low-energy (Kondo-like) physics is still domi- 
nated by the electron-electron interaction and room for 
nontrivial electron-phonon effects is left only for higher- 
frequency (atomic-like) processes. This has two very im- 
portant consequences: (i) the low-frequency slope of the 
dispersion is not strongly affected by the bare electron- 
phonon coupling and (ii) the deviation between the dis- 
persion with and without phonons (either in the form of 
a sharp kink or as a smooth change of slope) becomes 




FIG. 5: (Color online) PM calculation: (a) The Ek resolved 
spectral function for U = 6, n = 0.9 and A — 0.86. Note that 
for better visibility of all the peaks, the largest ones have been 
cut off at a certain value. This generates a small artificial 
plateau, otherwise absent, (b) The dispersion relation Ek as 
a function of £k for U — 6 and A = and A = 0.86. The 
horizontal dotted (dashed) line gives the energy wo (o>o), and 
the dashed line is for free electrons, Ek = £fc. (c) The real 
part of the self-energy E R (w) (with E R (0) subtracted) as a 
function of u for U — 6 and various values of A. (d) The 
effective mass as a function of A. 



a high-frequency electronic property, in striking contrast 
with the standard uncorrelated electron-phonon picture 
in which this is bound to the (renormalized) phonon fre- 
quency. 

The w-dependence of the real part of the self-energy, 
shown in panel (c), consequently shows relatively little 
modification at low energy and visible one at higher en- 
ergy. The behavior of the low energy slope and corre- 
spondingly effective mass is shown in part (d) as func- 
tion of A. As we noted above, there is an initial de- 
renormalization tendency upon increasing A, followed by 
a strong increase of m* for larger values of A. 

We noted for Fig. [S] (b) that no kink feature can easily 
be identified. This changes when we increase A. In Fig. 
[5] we show the dispersion at a larger value A=1.34. 
A kink feature is visible at ajfcsaO.15 while the renor- 
malized phonon frequency oj5~0.063 (extracted from the 
maximum of the phonon spectral function Pb(w) see in- 
set Fig. [3j. This clearly shows how these two quanti- 
ties, always bound to go hand in hand in a conventional 
uncorrelated metal as two aspects of the same physical 
phenomenon, can instead be different in a system which 
is described within paramagnetic DMFT. The renormal- 
ized phonon frequency wj still describes the influence of 
electron-phonon coupling on the phonon spectral func- 
tion, while the kink position lj^ is associated to the en- 
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FIG. 6: (Color online) PM calculation: The dispersion rela- 
tion Ek as a function of £*, for U = 6, n = 0.9 and different 
values of electron-phonon coupling. Included is a dashed line 
to fit the high energy slope of the dispersion. 



ergy range in the electronic spectrum where phonon ef- 
fects become visible. The latter is a new energy scale of 
the problem, which comes out of the competition between 
the instantaneous Hubbard repulsion and the retarded 
phonon- mediated attraction. It is interesting to see that, 
as seen in the inset of Fig. [31 pb(u>), which for A = is 
a delta function peaked at ujq, possesses a considerable 
amount of spectral weight for energies larger than Uq for 
this situation. By numerical fitting the magnitude of the 
kink was extracted to be m^ f=a 1.55. When we increase 
A further to overcome U, strong polaronic behavior will 
eventually dominate the low energy behavior again. 

Another possibility to find a kink due to electron- 
phonon coupling in the presence of large U is to increase 
the doping, as the effect of U to freeze the charge fluctua- 
tions is weakened then. We give an example for A = 0.86 
and 20% doping in Fig. [7J 

A weak kink feature can be identified at higher energy. 
This is more pronounced in calculations for larger values 
of A, but we stick to this value for comparison to later 
results. The magnitude of the kink is by numerical fitting 
rnfe w 1.39. Notice that also here the kink position is 
larger than the renormalized phonon energy. When we 
further increase the doping m^ will become larger and 
u>k will, for n ~ 0.5, eventually coincide 22 with uiq as 
expected in the conventional picture. 

We cannot compare the DMFT results obtained so far 
to experiments on cuprates as we have not yet intro- 
duced any effect coming from antiferromagnetic corre- 
lations. These are known to play a crucial role for small 
doping values, therefore we have to include them in our 
model calculation. One way of doing this is to let the 
system have long-range order so that we can describe (a 
classical mean-field version of) the physics of antiferro- 
magnetic exchange. This is what is done in the following 
section. 
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FIG. 7: (Color online) PM calculation: The dispersion rela- 
tion Ek as a function of ek for U — 6, n — 0.8 and different 
values of electron-phonon coupling. Included is a dashed line 
to fit the high energy slope of the dispersion. 



III. KINKS AND ANTIFERROMAGNETIC 
CORRELATIONS 

In order to include antiferromagnetic correlations we 
assume a bipartite lattice with A and B sublattice and, 
as in the standard DMFT implementation for long-range 
ordered phase a 18 ' 35 " — , we write the matrix Green's func- 
tion in the form 



G fc:CT M = 



1 



t-B^M £ k 



(a :<7 {u))( b ,o(oj) ~ £ 2 k \ £k UA"\ 



(4) 



with Cq,o-(w) = w+/x a](T — Y, airT (u>), a = A,B. For the AF 
order one has ha,o — A* — o"^s, ^B,a — M + crh s , and the 
condition Kg ]ff (w) = Y,a.-o{u). We consider solutions 
where the symmetry breaking field vanishes, h s — > 0. As 
we do not have any next-nearest neighbor hopping, the 
antiferromagnetic solution is lower in energy than the 
paramagnetic one for small values of the doping 35 . 

The form of the dispersion as found from the poles of 
the Green's function in Eq. @ can be obtained from the 
equation 



U± = -jU - 



E + ( w )± x /(A/2-S-H)2 + £ 2 



(5) 



where we only included the real parts of the self-energies. 
We defined fl = (fi-j- + Mi)/2, A/i = (/2-j- — fi,\)/2 with 
Mo- = Mo- - S(t,o, and we have split into static and dynamic 
part in the self-energy, ReE^u;) = E cr .o + E -(a;). We also 
introduced E ± (w) = (S t (w) ± E;(w))/2. The solution of 
Eq. |5]) as an implicit equation gives the dispersion Ek,± 
in the AF case. The mean field solution without phonons 
gives the two Slater bands, and upon hole doping the 
Fermi level moves into the lower one. By making a simple 
ansatz for the electron-phonon part of the self-energy, 
£9(w) = -\ a oj9(uj - |w|) with A = (A t + AJ/2, one 



can show that the different form of the equation does not 
lead to the occurrence of larger kinks in the dispersion 
relation than in the PM case, Eq. ([3]). 

Several studies pointed out the importance of coop- 
erative effects between lattice- and spin-polaron o 37 ' 45 , 
influencing both photoemissio n 28 i 29 i 46 ~ — and optical 
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of correlated systems. By studying an- 



tiferromagnetism in the Hubbard-Holstein model within 
DMFT it was clarified that a finite doping weakens these 
cooperative effects 23 . However, the values of A at which 
polarons are obtained in the antiferromagnetic phase are 
not only much smaller than the values needed to form 
local bipolarons but are also smaller than those needed 
to observe polaronic features in paramagnetic DMFT 
spectra^ 3 -. It is interesting to note that Dynamical Clus- 
ter Approximation (DCA) calculations 24 confirmed the 
picture found in AF-DMFT, even though short-range AF 
correlations in DCA are realistically present already at 
the paramagnetic level. 

The present AF-DMFT NRG study has the big advan- 
tage that we can resolve low-energy features in the spec- 
tral function at finite doping in a way that is beyond any 
of the previous exact diagonalization and Monte Carlo 
studies. In Fig. |8]we show the spin-resolved local spec- 
tral function in AF-DMFT solved with NRG. The pa- 
rameters, [7=6 and n=0.9, are the same as in Fig. @]and 
[5{b) where the spectral function without antiferromag- 
netic correlations was shown. As one can see the Fermi 
energy has moved into the lower band and at low energy 
the typical square root divergence in p^-(uj) appears. At 
the same time the minority spin pj,(w) displays the square 
root decrease. Hubbard band features at higher energies 
can also be identified. 

The sublattice spectral functions of the AF solution at 
10% doping and U = 6 show qualitatively similar features 
for A = and A = 0.86 as seen in Fig. [3] The spectrum 
at A=0.86 displays a strong quasiparticle excitation at 
the gap edge, with no big sign of renormalization. Yet 
the effect of the electron-phonon interaction is still much 
stronger than in the paramagnetic DMFT at the same U . 
To see this we show in Fig. [9] the dispersion relation for 
the same set of parameters, comparing again A=0 with 
A=0.86. 

A kink in E k is visible, while (see Fig. [SJd) just a smooth 
change in the slope could be observed for the same value 
of A without the AF correlations. By numerical fitting we 
can extract a magnitude of the kink of mfc«1.64. When 
doping is decreased or A increased the kink becomes even 
more pronounced. Similar values of rrik within the para- 
magnetic phase are reached only taking A=1.34. There- 
fore, the qualitative difference between paramagnetic and 
AF-DMFT coming from the underlying spin background 
reflects in a big quantitative difference in the strength of 
the kink. This is expected from the above arguments, but 
has not been shown explicitely for the dispersion relation. 

Particularly interesting for the present study is the 
kink position w^ in AF-DMFT. Contrary to the para- 
magnetic phase studied in Sec. Ill Bl the kink in the elec- 
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FIG. 8: (Color online) AF calculation: The local electronic 
spectral density for majority electrons pt(uj) and minority 
electrons p±(lj). 
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FIG. 9: (Color online) AF calculation: The dispersion relation 
Eh as a function of ek for U = 6, n = 0.9 and different values 
of electron-phonon coupling. Included is a dashed line to fit 
the high energy slope of the dispersion. 



tronic dispersion is at (or even slightly before) the renor- 
malised frequency u>q ~ 0.16, as shown in Fig. [9] This is 
very different from the paramagnetic DMFT result of the 
previous section, for which uiu looses its connection with 
the renormalized phonon frequency and becomes a higher 
frequency electronic property. For the values of A or dop- 
ing range where a clear kink can be identified (A=1.34, 
n < 0.8), indeed, paramagnetic DMFT gives a renormal- 
ized phonon frequency ujq < uju- Therefore we conclude 
that paramagnetic and antiferromagnetic DMFT calcu- 
lations can differ qualitatively in the physics behind the 
kink: in paramagnetic DMFT a larger A or larger dop- 
ing is needed to form a kink and the kink energy is not 
pinned to the renormalized phonon frequency, but rather 
is connected to a higher energy scale at the crossover 
between the Kondo-like physics and a more atomic-like 
one. In AF-DMFT a smaller value of A is enough to see 
a kink and the position of this kink closely tracks the 
renormalized phonon frequency. 

The stronger electron-phonon effects on the AF elec- 
tronic dispersion are connected directly to the real part 
of the self-energy, shown in Fig. [TU] for the majority (f) 
and minority (J.) spin component on the sublattice. 
Without the electron-phonon coupling minority electrons 
scatter much stronger with the majority spin electrons 
and therefore the slope of HP(ui) is larger than the one of 
£:?(w) for A = 0. 36 When A is increased, we can see that 
the low energy slope changes relatively little for S?(w). 
At A > 0.86 a feature at u ~ ujq can be identified in 
S^(w), which gives rise to the kink in the dispersion in 
Fig. O In calculations for larger A, smaller doping or 
smaller U this is more pronounced. This is related to 
the cooperative effect of the magnetic polaron interact- 
ing with phonons which was discussed earlier. The low 
energy slope of E?(ui) decreases on increasing A and no 
clear feature is seen around to ~ ujq . 

For small values of U and A these results for the self- 
energies can be understood from a perturbative calcula- 
tion to order g 2 (see Fig. Q}. As a starting point one can 
use the doped AF state and calculate the lowest order 
self-energy corrections due to coupling to the phonons. 
Then it is found that the majority spin electrons scatter 
much stronger with the phonons as the minority ones, 
and are therefore much stronger renormalized. H2(uS) 
shows a much more pronounced logarithmic divergence 
at Uq than in the normal state discussed in Sec. [TTJ For- 
mally, this comes from the enhanced density of states 
at the Fermi level which is a consequence of the AF or- 
der. An intuitive understanding can be motivated by 
considering that majority spin electrons on a sublattice 
are surrounded on average by the opposite spin direc- 
tion. This provides an energetic gain, but also leads to 
some degree of localization, as an itinerant excitation de- 
stroys the antiferromagnetic order. This effect can coop- 
erate with the coupling to the phonons inducing a fur- 
ther renormalization and as a result majority spins are 
more affected than without the ordered background. At 
the same time minority spin electrons couple less to the 
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FIG. 10: (Color online) AF calculation: The real part of the 
selfenergies £?"(w) (upper panel) and S^(oj) (lower panel) 
with Ect(0) subtracted. 



phonons. In addition to the different renormalization ef- 
fects, the electron-phonon coupling can lead to a stabili- 
sation and even an increase of the magnetization. These 
weak coupling considerations qualitatively carry over to 
the stronger coupling case relevant for the results above. 
However, additional renormalization effects due to U play 
a role then and make a perturbative treatment difficult. 



IV. CONCLUSIONS 

We have described the occurence of kinks in the elec- 
tronic dispersion relation due to coupling to a local 
phonon mode taking into account strong electronic cor- 
relations. As anticipated our results show that the sim- 
ple theory of kinks is modified for both strong electron- 
phonon and electron-electron coupling. For the paramag- 
netic state and small doping from half filling a clear kink 
is not visible anymore for large values of U>W unless 
the electron-phonon coupling assumes large values. The 
phonons can have an effect at higher energy whilst the 



low energy features are less affected. An important find- 
ing is that the effect of U to suppress the occurrence of 
kinks is less efficient if we allow for antiferromagnetic cor- 
relations. Then a cooperative effect leads to a stronger ef- 
fective coupling of the electronic quasiparticles (magnetic 
polarons) to the phonons. That something like this can 
happen has been observed in earlier studie o 23 ' . Here we 
clarify how this manifests itself for kinks. Nevertheless 
for the type of phonons we study in competition with 
the Coulomb interaction we find that inspite of relatively 
large values, A ~ 0.8, kink features are moderate even 
with antiferromagnetic correlations, rrik — 1-6. This may 
be a consequence of the type of phonons we study (Hol- 
stein) . Some works have indeed shown that non-Holstein 
phonons have stronger polaronic effects in correlated sys- 
tems, due to a less direct competition with the local Hub- 
bard repulsion^—. Therefore a quantitative comparison 
with results for the low-energy kink in cuprate supercon- 
ductors, or a statement about the origin of the nodal kink 
in these systems is hard to make here. 

Still, we can relate our findings to the puzzling results 
for the kink position reported by Kordyuk et al. in Ref. 
[13 . There, as already discussed in the Introduction, a 
suprisingly strong dependence on doping and tempera- 
ture of the kink position was observed. Nothing like this 
is can be obtained within conventional electron-phonon 
theories. Our results show that, by taking the effect of 
electronic correlation properly into account, the kink po- 
sition becomes strongly dependent on the strength of an- 
tiferromagnetic correlations. We found that when AF 
correlations are strong, the kink is close to the renor- 
malized phonon frequency (i.e. is more "low energy"), 
while when AF correlations are absent, the kink can oc- 
cur at higher frequencies and also needs a larger electron- 
phonon coupling to be of similar strength. Therefore it 
is very plausible to expect that, at low doping, where 
AF correlations are stronger, the kink energy is small in 
absolute value (see Fig. [5]) , while at larger values of the 
doping, where AF correlations are weaker, the kink oc- 
curs at larger energies (see Fig. [7]). Interestingly, this is 
exactly the trend shown in the data on the Bi-compound 
at the lowest temperature^. 

Generally, it is very important to test - as we do here - 
the predictions of theories which, like paramagnetic and 
AF-DMFT, describe electron-phonon and strong elec- 
tronic correlations on an equal footing. Strong coupling 
to a phonon mode is very often discarded as being re- 
sponsible for the nodal kink in cuprates on the basis of 
cliches coming from conventional theories neglecting elec- 
tronic correlations. The doping dependence is one typical 
example. Our study shows instead that, when including 
correlation, properties like the doping dependence of the 
kink energy may become highly non-trivial, and there- 
fore the conventional picture cannot be used to rule out 
phonons. 
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